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1 .  Summary  of  Solar  Seismology  Work  Performed 
During  1982-1986  and  Supported  by  APOSR 
Contract  82-0019 


1.1  Scope  of  the  1982-1986  Research  Program 

In  the  1970s  an  observational  technique  for  the  study  of  solar 
oscillations  with  periods  longer  than  5  min  was  developed  at  SCLERA.  The 
technique  yields  information  about  the  temperature  eigenfunction  of  an 
oscillation  at  the  extreme  solar  limb.  The  primary  goal  of  the 
observational  program  for  1982-1986  was  to  use  this  technique  to 
systematically  increase  the  data  base  for  the  study  of  solar  oscillations. 
Because  of  the  complexity  of  the  eigenfrequency  spectrum  and  because  of  the 
relatively  small  signals  produced  by  the  oscillations,  the  analysis  of  the 
observations  for  identification  and  classification  of  modes  of  oscillations 
has  proved  to  be  very  difficult.  The  main  objective  of  the  data  analysis 
program  for  1982-1986  was  to  develop  analysis  techniques  which  would  lead  to 
mode  classifications  based  on  the  SCLERA  differential  radius  observations 
and  on  the  combination  of  the  differential  radius  observations  with 
differential  velocity  observations.  The  goal  of  the  theoretical  program  has 
been  the  inversion  of  the  multiplet  fine  structure  obtained  in  the  data 
analysis  program  in  order  first  to  test  it  for  internal  consistency  and 
second  to  infer  information  about  the  radial  and  latitudinal  dependence  of 
the  Internal  rotation  of  the  Sun. 

The  results  of  these  programs  appear  in  a  series  of  works  and  have  led 
to  the  publication  of  a  sizable  amount  of  data;  the  development  of  the 
analysis  programs  for  the  classification  of  modes;  the  classification  of 
•  1209  resolved  modes  belonging  to  166  multlplets,  including  acoustic,  f-, 
and  gravity  modes;  tests  of  mode  classification  programs;  Inferred 
information  on  the  internal  rotation  of  the  Sun;  and  evidence  of  mode 
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coupling  which  nay  be  relevant  to  the  solar  neutrino  paradox.  These  results 
are  reviewed  in  the  following  sections. 

1.2  Observations  and  Publication  Schedule  of  Observational  Data 

1.2.1  Schedule  of  Observations 

The  1983  and  1 98**  Observations 

A  series  of  observations  were  made  in  1983  and  1984  for  which  the 
observational  technique  and  telescope  apparatus  were  very  similar  to  those 
used  in  the  1979  SCLERA  differential  radius  observations  (cf.  Bos  and  Hill 
1983)*  In  this  work,  data  was  obtained  for  the  solar  latitudes  -45°.  0s, 
45*.  and  90*.  Data  was  acquired  at  64-second  intervals  for  each  of  these 
four  latitudes,  and  the  entire  256-second  cycle  was  repeated  continuously 
throughout  each  observing  day.  The  set  of  data  was  obtained  between  May 
1983  and  March  1984.  From  6/1/83  to  7/20/83  data  was  taken  on  47  days  at  an 
average  of  8.8  hours  per  day.  From  10/10/83  to  3/21/84  data  was  taken  on  29 
days  at  an  average  of  5.8  hours  per  day. 

The  1985  Observations 


The  observations  made  in  1985  were  obtained  with  a  new  observational 
technique  described  in  Section  1.7.  With  this  technique,  solar  oscillations 
are  detected  as  perturbations  In  the  radiation  intensity  at  disk  center  as  a 
function  of  x  for  0.5  <  X  <  1.7  u.  Observations  were  made  on  35  days 
spanning  a  70-day  period  from  March  23  to  May  31 .  1985  at  an  average  of  8.3 
hours  per  day. 

The  1986  Observations 


The  observational  technique  and  telescope  apparatus  for  the  1986 
observations  were  also  very  similar  to  those  used  in  the  1979  dlf'erential 


2 


radius  observations  (cf.  Bos  and  Hill  1983)  with  emphasis  placed  on 
equatorial  observations.  The  telescope  was  continuously  manned  during  the 
spring  and  early  summer  observing  season  of  1986.  Although  the  particular 
observing  season  proved  to  be  exceedingly  cloudy,  over  100  hours  of  data 
were  obtained. 

1.2.2  Publication  Schedule  of  SCLERA  Observational  Oata 

Data  from  SCLERA  observations  have  been  presented  in  a  series  of  works 
for  the  1978  and  1979  observations.  The  average  of  13  dally  power  spectra 
with  a  30  uHz  resolution  is  shown  in  Figure  1  of  Caudell  et  al.  (1980)  and 
also  in  Figure  10  of  Caudell  (1980)  for  the  1978  observations.  Also,  for 
the  1978  observations,  the  diagrams  of  phase  versus  time  for  12  peaks  in  the 
power  spectrum  are  shown  in  Figures  3,  4,  and  5  of  Caudell  et  al.  (1980)  and 
also  in  Figures  15,  16,  and  17  of  Caudell  (1980).  For  the  1979 
observations,  the  list  of  works  containing  data  is  more  extensive  and  is 
given  in  Table  1.  The  data  in  these  works,  much  of  it  available  as  early  as 
1982,  permits  certain  of  the ’results  to  be  Independently  tested  such  as  was 
done  for  the  1978  observations.  In  that  test,  Gough  (1980)  Independently 
examined  the  phase  diagrams  for  the  1978  observations  and  concurred  in  the 
conclusions  of  Caudell  et  al.  (1980).  The  first  data  from  the  1983  and  1985 
observations  to  be  published  are  scheduled  to  be  available  in  1986  (Yi  1986 
and  Oglesby  1986,  respectively). 
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Table  1 

SCLERA  Publications  of  Data  to  Date  for  the  1979  Observations 


Source 

Fig.  Number 

Resolution 

Type  of  Power 

Frequency  Range 

Spectra® 

Bos  (1982) 

4.1 

30  uHz 

0  -  3.2  mHz 

4.2 

vt 

ft 

“.3 

ft 

P4 

tt 

4.4 

It 

pJ 

ft 

3 

4.5 

0.28  uHz 

;i 

240  -  272  uHz 

4.6 

t* 

ft 

2 

4.7 

0.28  uHz 

V 

r  o 

420  -  452  uHz 

4.8 

H 

tt 

4.9 

H 

tt 

4.10 

II 

tt 

4.11 

0.28  kHz 

Pp. 

*2 

550  *  582  uHz 

4.12 

i« 

If 

0.28  uHz 

P1 

420  -  452  uHz 

6.3 

0.28  uHz 

P1 

235  -  285  uHz 

Hill,  Bos, 

&  Goode 

(1982) 

1 

0.28  uHz 

P1 

235  -  285  uHz 

Bos  &  Hill 

(1983) 

3a 

0.28  yHz 

P1 

450  -  482  uHz 

Hill  (1984a) 

1b 

0.28  uHz 

P1 

235  -  285  uHz 

Hill  (1985) 

1 

0.28  uHz 

pP1 

2492  -  2502  uHz 

2 

ft 

3227  -  3237  uHz 

3 

tl 

P2 

3909  -  3919  uHz 

4 

ft 

?2 

3209  -  3225  uHz 

5 

If 

?2 

3616  -  3632  uHz 

6 

ft 

P2 

3890  -  3906  uHz 

Hill,  Taah,  & 
Pad in  (1986) 

1 

0.28  uHz 

P1 

100  -  110  uHz 

Hill  & 

Czar nows ki 

3 

0.28  uHz 

P1 

r2 

0  -  20  uHz 

(1986) 

4 

If 

ft 

£saae  as  Fig. 

6.3  of  Bos  >  • 

??2)  and  Fig.  1 

of  Hill  (1984a). 

Saae  as  Fig. 

6.3  of  Bos  ( ' 

and  Fig.  1 

of  Hill,  Bos,  and 

Goode  (1982). 

Definition  of 

P,  noaenclat 

.-e  «ane  as  used  by  Hill  (1984b). 
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1.3  The  SCLERA  Mode  Classification  Program 

The  objective  of  the  mode  classification  program  developed  at  SCLERA  is 
the  classification  of  low  to  intermediate  degree  acoustic  and  gravity  modes 
with  low  to  intermediate  radial  orders.  This  program  is  based  primarily  on 
the  differential  radius  observations  from  SCLERA  supplemented  with  various 
types  of  Doppler  shift  and  total  irradiance  observations  from  other 
observatories .  These  observations  provide  the  following  categories  of 
information  for  use  in  mode  classification:  (1)  the  multiplet  fine  structure 
of  the  eigenfrequency  spectrum  expected  for  a  slowly  rotating  axisymmetrlc 
system,  (2)  the  symmetry  properties  of  the  eigenfunction,  (3)  the  parities 
of  i  and  m,  (4)  the  magnitude  of  1,  (5)  the  (exp  lme)~dependence  of  the 
eigenfunction,  and  (6)  the  e-dependence  of  the  eigenfunction  given  by  Y®. 
In  addition,  a  theoretical  eigenfrequency  spectrum  from  a  standard  solar 
model  is  used  in  the  classification  of  the  radial  order. 

The  subsections  below  review  how  the  above  information  is  obtained 
and/or  applied  to  mode  classification.  It  is  apparent  that  this  rather 
extensive  list  may  permit  not  only  the  classification  of  modes  but  also  the 
implementation  of  a  number  of  tests  of  the  accuracy  of  a  given  set  of  mode 
classifications.  The  availability  of  a  set  of  independent  tests  has  become 
a  hallmark  of  the  SCLERA  mode  classification  program.  These  tests  are  in 
addition  to  those  described  in  Section  1.5  concerned  with  the 
reproducibility  of  results  from  one  year  to  another. 

1.3.1  Symmetry  Properties  of  Eigenfunctions 

The  determination  of  the  symmetry  properties  of  the  eigenfunctions  has 
been  one  of  the  most  important  exercises  of  the  last  few  years.  First,  if 
the  Sun  were  not  an  axisymmetrlc  system  as  seen  by  the  normal  modes,  the 
complexity  of  the  multiplet  fine  structure  could  considerably  complicate  the 
mode  classifications.  Thus  the  information  obtained  on  the  eigenfunction 
symmetry  properties  was  first  used  to  determine  the  type  of  multiplet  fine 
structure  to  be  expected  (see  Bos  and  Hill  1983)  • 

With  the  determination  of  the  symmetry  properties  of  the  system,  two 
analysis  techniques  based  on  the  inferred  symmetry  properties  of  the 
eigenfunctions  have  been  successfully  implemented.  The  first  technique  uses 
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the  symmetry  properties  to  generate  power  spectra  which  contain  essentially 
the  signals  for  modes  with  only  one  combination  of  tne  parities  of  i  and  m. 
The  four  different  combinations  of  the  parities  of  i,  and  m,  respectively, 
are  (even,  even),  (even,  odd),  (odd,  even),  and  (odd,  odd).  This  technique 
leads  to  an  important  reduction  in  the  density  of  peaks  in  the  complex  power 
spectra  and  permits  the  immediate  classification  of  the  associated  parities 
of  l  and  mwith  a  relatively  high  probability  of  being  correct.  Examples  of 
the  successful  use  of  this  technique  are  found  in  the  works  of  Bos  and  Hill 
(1983).  Hill,  Bos,  and  Goode  (1982),  Hill  (1984b),  and  Hill  (1985a, b). 

The  second  technique  uses  the  symmetry  properties  to  test  for  the 
cocrectness  of  a  particular  mode  classification.  If  the  mode 
classifications  for  multiplets  are  correct,  the  signals  in  the  complex  power 
spectra  due  to  every  other  member  of  a  multiplet  should  be  confined  to  one 
power  spectrum  while  the  signals  due  to  the  remaining  set  should  be  confined 
to  a  second  power  spectrum.  Examples  of  this  technique  are  found  in  the 
same  works  as  listed  in  the  previous  paragraph.  It  should  be  noted  that  the 
second  technique  permits  the  determination  of  what  fraction  of  the  time 
Incorrect  peak  classifications  are  made  due  to  coincidental  alignment  of  the 
frequency  of  one  mode  with  that  of  an  alias  of  a  second  mode  (cf.  Section 
I V£  of  Hill  1984b  and  Section  IVe  of  Hill  1985b). 

1.3*2  Frequency  Patterns  of  Multiplet  Fine  Structure 

The  determination  of  the  symmetry  properties  of  the  eigenfunctions  as 
those  arising  from  an  axisymmetric  system  (Bos  and  Hill  1983,  Hill  1984b) 
furnished  an  extremely  Important  foundation  for  defining  the  basic 
properties  of  the  multiplet  fine  structure.  In  particular,  a  slowly 
rotating  axisymmetric  Sun  should  lead  to  a  multiplet  fine  structure  that  to 
a  first  approximation  is  uniformly  spaced  in  frequency.  This  frequency 
pattern  is  often  referred  to  as  a  Zeeman  pattern,  and  its  simplicity  has 
played  an  important  role  in  the  SCLERA  mode  classification  program. 
Examples  of  the  value  of  this  technique  are  found  in  the  first  five 
subsections  of  Section  1.4. 
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1.3*3  Azimuthal  Dependence  of  Eigenfunctions 


The  azimuthal  dependence  (exp  im$)  of  an  eigenfunction  expected  for  an 
axisymmetric  system  is  the  basis  of  an  important  technique  developed  by  Hill 
(1984b)  and  extended  in  the  work  of  Hill,  Alexander,  and  Caudell  (1985). 
The  information  on  the  m-dependence  is  contained  in  the  observable  D^j  , 
which  is  a  functional  of  (exp  im*).  Formally,  is  essentially  the  ratio 
of  the  finite  Fourier  transform  of  the  sin  m*  component  of  (exp  im+)  to  the 
finite  Fourier  transform  of  the  cos  m+  component  of  (exp  im$) .  The  general 
properties  of  are  that  it  is  a  quasl-periodic  function  of  m  and  is 
antisymmetric  for  m-*-m.  An  example  of  the  is  shown  in  Figure  1,  which 
was  obtained  in  the  work  on  the  low-order,  low-degree  acoustic  modes  (cf. 
Section  1.4.2). 

The  functional  has  been  important  in  the  determination  of  the  m  -  0 
member  of  a  multiplet  because  of  its  symmetry  properties  about  m  •  0.  It 
has  also  been  an  extremely  valuable  part  of  the  SCLERA  mode  classification 
program  because  it  furnishes  a  highly  statistically  significant  test  of  a 
particular  set  of  mode  classifications.  Examples  of  its  application  are 
found  in  the  first  five  subsections  of  Section  1.4. 

1.3.4  Zeros  of  Y® 

A  second  technique  has  been  introduced  to  determine  the  m  -  0  mode  and 
the  A  value  of  a  multiplet  for  modes  with  |m|  >  20.  This  technique  is  based 
on  the  properties  of  Y^  where  is  the  average  of  Y®  over  the  detector 
geometry  used  by  Bos  (1982)  (see  Bos  and  Hill  1983).  The  relationship 
between  l  and  the  values  of  m,  m^,  for  which  Y®  •  0,  are  shown  in  Table  2 
for  even  intermediate  values  of  l. 

The  first  successful  use  of  these  properties  with  respect  to  zeros  was 
in  the  work  on  phase-locked  modes  (cf.  Section  1.4.5).  For  such  modes,  the 

zeros  of  Y^  introduce  1 80 0  phase  shifts  in  the  plots  of  phases  of  the  modes 
versus  m. 
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Fig  1  The  e  *  dependence  of  the  eigenfunctions  as  represented  by  <0^. 
The  function  D  is  defined  in  terms  of  the  ratio  of  the  finite 

*  J 

Fourier  transforms  of  the  sin  m*  and  cos  m+  dependent  terms  of  the 
oscillatory  component  of  the  radiative  intensity'^'.  The  modes 
Included  in  this  analysis  by  Hill  ( 1 984b )  are  the  p1 

P23 . P2,10:  and  p3,2’  p3, 3’  p3, m*  The  solid  curve  18  the 

predicted  <D,.>  using  the  inferred  I'  obtained  by  Hill,  Alexander, 
1 J  ’  o 

and  Caudell  (1985).  The  dashed  curve  is  the  predicted  <Dij>  for  a 
theoretically  derived  1'. 
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Table  2 


Values  of  d,  mQ,  for  y“  -  0 


l 

mo 

26 

9 

28 

13 

30 

17 

32 

21 

31* 

25 

1.3.5  The  Ratio  of  Aaplltudes  from  Differential  Velocity  and  Differential 
Radius  Observations 

The  gravity  mode  spectrum  is  much  more  complex  than  the  low-order,  low- 
degree  acoustic  mode  spectrum.  Evidence  for  the  detection  of  gravity  modes 
was  presented  by  Hill  and  Caudell  0  979)  and  by  Bos  and  Hill  0983),  but  the 
classification  of  the  modes  has  been  hampered  by  the  complexity  of  the 
gravity  mode  spectrum.  Hill  and  Caudell  0  979)  were  able  to  estimate  the 
magnitude  of  l  only.  Hill,  Bos,  and  Goode  0982)  identified  a  subset  of  two 
gravity  mode  multiplets,  but  were  only  able  to  place  a  lower  limit  on  l 
using  Just  the  differential  radius  observations.  However,  the  complexity 
can  be  managed  with  the  additional  information  gained  on  the  degree  of  a 
mode  as  a  result  of  the  positive  results  obtained  by  Hill  (1985b)  in  a  test 
of  the  Hill,  Tash,  and  Padin  0986)  apparent- velocity  hypothesis.  In  this 
circumstance,  the  observed  constraints  on  l  and  m  are  decoupled. 

The  development  of  the  apparent-velocity  hypothesis  described  in 
Section  1.6.2  arose  out  of  a  comparison  of  results  obtained  with  different 
observing  techniques.  This  hypothesis  interprets  the  velocity  signal 
obtained  in  differential  Doppler  shift  studies  of  long-period  oscillations 
as  an  apparent  velocity  Cue  to  the  combined  effect  of  the  surface  rotation 
of  the  Sun  and  a  pert ur sat  ion  in  the  radiation  intensity  produced  by  a 
normal  mode  of  oscillation.  The  ratio  of  this  velocity  signal  V  to  the 
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observed  amplitude  of  differential  radius  amplitude  ar1  can  be  written 
according  to  Hill,  Tash,  and  Padin  (1986)  as 


mar 

vC— 7— 5-3 

w(m) 


/ir/2f[g» 
J0  L80 


-(i)“a. 
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Y®(tt/2,  0) 

‘a  KH)  *  ~ 

YlU/2,  0) 
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1.067g2  Sin2  0^  sin2  0  003  0  d0. 


Where  ag  is  a  proportionality  constant  for  the  symmetric  line  contribution, 
Y®(*/2,0)  is  the  average  of  y”  over  the  detector  geometry  used  by  Bos  and 
Hill  (1983),  g^  and  g£  are  the  effective  first  and  third  derivatives  with 
respect  to  wavelength  of  $  '  in  the  spectral  line,  f  is  the  spatial  filter 
function  of  the  velocity  detector,  and  w(m)  is  the  spatial  filter  function 
for  the  differential  radius  observations  given  by  Hill,  Alexander,  and 
Caudell  (1985). 

The  ratio  of  m/w(m)  is  only  weakly  dependent  on  m  for  2  <  |m|  <  10  so 
that  the  ratio  V/Cmar^/wCm]  is  to  a  good  approximation  only  a  function  of  l. 
Examples  of  what  might  be  expected  are  shown  in  Figure  2.  Thus,  Information 
is  available  on  the  magnitude  of  l  by  combining  the  results  of  differential 
velocity  and  differential  radius  observations.  An  example  of  the 
application  of  this  to  mode  classification  is  found  in  the  work  <*f  Hill 
(1985b)  described  in  Section  I.4.M. 


1.3.6  The  Ratio  of  Amplitudes  from  Total  Irradiance  and  Differential  Radius 
Observations 


A  technique  has  been  developed  to  obtain  detailed  information  on  the 
1  lmb-dar kerning  properties  of  the  oscillatory  component  of  the  total 
radiation  intensity  associated  with  an  oscillation.  This  amplitude  of  the 
total  radiation  intensity  can  be  written  as 

J'  -  i» ( w)y“( e, ♦)  (2) 

where  u  is  the  cosine  of  the  heliocentric  angle  between  the  observer  and  the 
direction  defined  by  (0,$)  and  I'(u)  is  only  weakly  dependent  on  l.  In 
particular,  it  is  the  limb-darkening  properties  of  l'(u)  that  are  of 
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Fig  2  The  predicted  ratio  of  the  symmetric  line  apparent  velocity  signal 

to  the  differential  radius  signal  for  different  velocity  detector 

configurations.  The  predicted  apparent  velocity  is  based  on  the 

properties  of  the  Fel  51 2M  transition.  Two  curves  are  for  the 

differential  velocity  detector  configurations  used  at  the  Crimean 

and  Stanford  observatories .  The  predicted  apparent  velocity  for 

the  whole  disk  observations  assumed  spectrometer  specifications 

which  are  the  same  as  those  of  the  Stanford  differential  velocity 

observations.  The  remaining  curve  is  for  a  disk  annulus 

configuration  used  in  the  Crimean  observations ,  but  with 

spectrometer  slits  centered  at  ±  0.035  A  off  line  center  with  slit 
0 

widths  of  0.030  A. 


11 


interest  here.  With  Equation  (2),  the  ratio  of  the  fractional  change  in 
total  irradiance  S£/Sp  to  the  differential  radius  amplitude  can  be 
written  as 


Ar1 

Csf/sf3CS7I)3 


/;  p”  gp  sinecosede 
gp  sinecosede 


(3) 


where  I'(u)  has  been  expressed  as 

I'<u>  -  I£g*.  (u>  (4) 

(see  Hill  and  Kroll  1986).  The  factors  Q  and  1^  are  constants.  The 
parameter  w(m)  is  the  factor  that  relates  the  observed  signal  ar^  to  the 
fractional  change  in  the  radiation  intensity  at  the  extreme  equatorial  limb 
and  is  given  by  Hill,  Alexander,  and  Caudell  (1985). 

The  information  on  the  spatial  properties  of  g'  is  contained  in  the 
Fredholm  Integral  equations  of  the  first  kind  given  by  Equation  (3).  This 
information  is  obtained  by  the  inversion  of  integral  equations  for  a  given 
set  of  i. 

This  technique  has  been  successfully  employed  by  Hill  and  Kroll  (1986) 

in  the  inversion  of  a  set  of  six  integral  equations  for  l  -  1 . 5.  The 

resulting  g£  is  shown  in  Figure  3. 

1.4  Summary  of  Modes  Classified  and  Observationally  Derived  Elgenfrequency 
Diagram 

2 

1.4.1  Low-Degree  5  min  Oscillations 

Individually  resolved  low-degree  5  min  modes  have  been  identified  using 
the  differential  radius  observations  from  SCLERA  augmented  with  the 
information  on  the  elgenfrequency  spectrum  obtained  from  Doppler  shift  and 
total  irradiance  observations  (Claverie  e£  al.  1981;  Grec,  Fossat,  and 
Pomerantz  1983;  Duvall  and  Harvey  1983;  and  Woodard  and  Hudson  1983).  In 


2.  This  is  a  synopsis  of  Hill  (1985a). 
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Fig  3 


The  radial  dependence  of  g£  obtained  by  Hill  and  Kroll  (1936)  in 
the  inversion  of  integral  Equations  of  the  type  given  by  Equation 
3.  The  observations  used  were  the  differential  radius 
observations  of  1979,  the  total  irradiance  observations  of  Woodard 
and  Hudson  (1983),  and  the  mode  classifications  by  Hill  (1985b). 
The  value  of  t  indicated  in  this  analysis  are  i.  -  1 . 5. 
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particular,  independent  works  on  the  5  min  oscillations  were  used  to  furnish 
measures  of  the  m  -  0  eigenfrequencies  for  the  low-degree  5  min  modes  as  a 
starting  point  in  classifying  the  peaks  in  the  power  spectra  of  the  1979 
differential  radius  observations.  With  this  information,  the  l  -  0  and 
t  -  1,  m  -  0  singlets  were  identified  first,  followed  by  the  1  •  1 ,  m  »  ±1 
doublets.-  Once  the  multlplet  fine  structure  was  obtained  for  the  1  -  1 
multi plet,  the  analysis  was  extended  to  include  the  l  •  2  multiplet.  This 
procedure  was  repeated  until  multiplets  with  1  »  6  had  been  classified. 

A  total  of  184  modes  were  classified  belonging  to  83  multiplets.  The 
multiplets  are  indicated  in  Figure  6  (cf.  Section  1.4.6).  It  should  be 
noted  .that,  in  this  analysis,  consideration  was  given  to  the  complexities 
Introduced  in  the  power  spectra  when  the  mode  coherence  time  is  less  than 
the  length  of  the  observing  period.  Also,  the  symmetry  properties  of  the 
signals  associated  with  the  classified  modes  furnished  a  statistically 
significant  teat  of  the  classifications. 

1.4.2  Low-order,  Low-degree  Acoustic  Modes ^ 

The  first  major  mode  classification  project  at  SCLERA  was  concerned 
with  the  low-degree  acoustic  modes  for  frequencies  Just  above  the  asymptotic 
gravity  mode  limit  (see  Figure  6  of  Section  1.4.6).  The  reduced  complexity 
of  the  elgenfrequency  diagram  in  this  region  and  the  intrinsic  line  widths 
for  these  modes  of  <  0.1  yHz  were  the  two  factors  that  dictated  this 
frequency  region  for  study  in  the  first  analysis.  It  was  in  this  project 
that  the  fundamental  features  of  the  SCLERA  mode  classification  program  were 
developed,  as  outlined  in  Section  1.3. 

Identification  of  multiplets  was  implemented  as  an  Iterative  process. 
Initially,  the  Zeeman-like  frequency  pattern  was  sought  for  a  multlplet 
where  the  pattern  was  most  easily  discernible  at  statistically  significant 
levels  and  the  associated  3  (see  Section  1.3.3)  exhibited  physically 
plausible  properties.  In  tne  multlplet  subsequently  classified  as  p1  2Q,  a 
relatively  large  fraction  of  tne  maximum  number  of  peaks  possible,  21  ♦  1  , 
was  tentatively  identified.  The  measures  of  the  were  used  for 


3.  This  is  a  synopsis  of  Hill  (1984). 
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Identifying  the  m  •  0  member  of  the  multiplet.  With  the  identification  of 
the  m  •  0  mode,  the  classification  of  n  was  made  by  comparing  the  observed 
m  ■  0  eigenfrequency  to  that  predicted  theoretically  with  the  standard  solar 
model  of  Salo  (1982).  When  the  most  likely  location  of  the  m  -  0  mode  was 
defined  and  the  multiplet  classified  according  to  n,  the  frequency  pattern 
of  the  multiplet  was  fit,  using  a  least-squares  analysis,  to  a  power  series 
cubic  in  m  written  as 


nlm 


v_.  ♦  vl.m  ♦  v?.m2/2  ♦  v"!m3/3! 


nl 


nl 


nl 


nt 


(5) 


where  v  .  is  the  eigenfrequency  of  the  mode  specified  by  (n,t,m).  The 
nun 

actual  final  selection  of  members  of  the  multiplet  was  based  upon  the 
deviation,  &v,  between  the  frequency  of  a  peak  under  consideration  and  the 
frequency  determined  by  the  polynomial  fit.  The  maximum  deviation  accepted 
was  ±0.06  yHz,  approximately  twice  the  standard  deviation  of  &v  (see  Section 
IVc  of  Hill  198Mb) . 

After  the  analysis  of  the  above  multiplet  was  completed,  identification 
of  contiguous  multiplets  in  l  was  undertaken,  information  about  the 
rotational  splitting  effects  linear,  quadratic,  and  cubic  in  m  and  about  the 
m  •  0  mode  eigenfrequency  obtained  from  the  least- squares  analysis  of  the 
multiplet  fine  structure  was  used  in  conjunction  with  the  theoretical 
eigenfrequency  spectrum  of  Saio's  (1982)  model  to  predict  the  parameters  of 
multiplets  contiguous  in  i.  With  the  inclusion  of  the  predicted  values  of 
these  parameters,  the  process  for  the  previous  multiplet  was  repeated, 
adding  one  constraint:  the  composite  list  of  observed  multiplet  parameters 
must  be  Internally  consistent. 

The  process  was  reiterated  until  candidates  for  n  -  1  multiplets  were 
identified  for  t  •  4,..., 22,  for  n  -  2  multiplets  for  t  •  3,..., 10,  and  for 
n  •  3  multiplets  for  t  •  2,  3,  and  4.  a  total  of  293  modes  were  classified 
belonging  to  30  multiplets.  These  multiplets  are  included  in  Figure  6  of 
Section  1.4.6. 
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1.4.3  Intermediate-Degree  f -modes 

The  mode  classification  project  for  the  low-order,  low-degree  acoustic 
modes  took  advantage  of  the  minimum  in  the  complexity  of  the  eigenfrequency 
spectrum  found  Just  above  the  asymptotic  g-mode  limit  (see  Section  1.4.2). 
For  this  same  frequency  region,  the  intermediate-degree  f-modes  (n  -  0)  are 
expected  based  on  the  theoretical  eigenfrequency  spectrum  of  a  standard 
solar  model.  These  modes  are  of  particular  Interest  because  they  can  play  an 
Important  role  in  studying  the  outer  -  10%,  by  radius,  of  the  solar 
Interior. 

Preliminary  evidence  of  intermediate-degree  f-modes  was  found  in  a  test 
of  the  validity  of  a  uniformly  rotating  convection  zone  model  by  Hill, 
Rosenwald,  and  Rabaey  (1985).  In  the  subsequent  work  of  Hill  (1986b).  which 
extended  the  analysis  of  Hill,  Rosenwald,  and  Rabaey  (1985),  the  multiplet 
fine  structure  for  the  intermediate-degree  f-mode  multi plets  of  angular 
degree  t  •  33  was  Identified.  With  this  finding  as  a  starting  point  in  a 
mode  classification  program  similar  to  that  used  by  Hill  (1984b),  the  study 
of  the  intermediate-degree  f-modes  in  the  450-650  uHz  region  of  the  spectrum 
was  initiated. 

In  the  classification  of  the  l  -  33  multiplet,  the  frequencies  of  modes 
were  fit  to  Equation  (5).  The  final  selection  of  members  of  the  multiplet 
was  based  upon  the  deviation,  Av,  between  the  frequency  of  a  peak  under 
consideration  and  the  frequency  determined  by  the  polynomial  fit.  After  the 
analysis  of  this  multiplet  was  completed,  the  analysis  of  contiguous 
multiplets  in  l  was  undertaken.  Information  about  the  rotational  splitting 
effects  linear,  quadratic,  and  cubic  in  m  and  about  the  m  -  0  mode 
eigenfrequency  obtained  from  the  least-squares  analysis  of  the  multiplet 
frequency  pattern  was  used  in  conjunction  with  the  theoretical 
eigenfrequency  spectrum  of  Saio's  (1982)  model  to  predict  the  parameters  of 
the  multiplets. contiguous  in  1.  With  the  Inclusion  of  these  parameters,  the 
process  for  the  first  multiplet  was  repeated  with  one  added  constraint:  the 
composite  list  of  observed  multiplets  must  be  internally  consistent  as 
measured  by  the  following  criteria.  These  criteria  are:  (1)  the  modes 


4.  This  is  a  synopsis  of  the  work  by  Rabaey  and  Hill  (1986) 
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exhibit  an  almost  uniformly  spaced.  Zeeman-like  frequency  pattern;  (2)  the 
modes  possess  minimum  quadratic  and  cubic  terms  in  the  departure  from  the 
uniform  spacing;  (3)  the  linear,  quadratic  and  cubic  terms  are  consistent 
with  those  fornd  for  neighboring  multiplets  in  t;  (4)  the  frequency  of  the 
axlsymmetric  mode  is  consistent  with  the  frequencies  of  corresponding  modes 
in  neighboring  multiplets;  and  (5)  the  classification  of  multiplets 
identified  by  criteria  (l)-(4)  is  confirmed  by  the  observed  horizontal 
spatial  properties  of  the  oscillations  as  measured  by  D^.  As  noted  in 
Section  1.3.3S  is  a  f motion  of  m  which  allows  the  confirmation  of  the 
m  -  o  mode  of  the  multiplet. 

With  this  process  a  set  of  18  multiplets  and  a  total  of  nearly  340 
modes  were  identified  and  classified  with  i  values  from  19  to  36  (Rabaey  and 
Hill  1986.)  The  resulting  m  «  0  elgenfrequencles  vQ  are  indicated  in 
Figure  6  of  Section  1.4.6.  In  Figure  4  the  results  of  the  linear  term  in  m, 
v*  ,  which  is  due  to  rotational  splitting  effects,  are  shown.  The  error 

0  •  It 

bars  shown  are  formal  errors  based  upon  the  least-squares  analysis.  Since 

the  rotational  splitting  kernels  for  the  intermediate-degree  f-modes  are 

well  localized  in  the  outer  101  of  the  Sun,  one  can  easily  verify  that  these 

values  of  v'  #  are  consistent  with  the  rotational  curve  obtained  by  Hill, 

Rosenwald,  and  Rabaey  (1985).  In  Figure  5  a  comparison  is  made  of  NQ 

and  p  where  NQ  ^  is  the  number  of  modes  identified  for  a  given  l,  Nmj{  Is 

the  maximum  number  expected  for  that  multiplet  and  p  is  the  probability  of 

finding  one  or  more  peaks  from  a  random  distribution  within  *0.065  uHz  of 

v_  .  obtained  in  the  polynomial  fits  (see  Hill  1984b).  The  values  of  p 

u,  L,m 

range  from  0.301  to  0.346  and  the  values  of  N-./N..,  range  from  a  low  of 
—  ui  max 

0.321  to  a  maximum  of  0.469.  Computing  d,  the  weighted  average  of  the 

difference  between  the  observed  ratio  N.  ,/N  _  and  p,  we  find  d  •  0.088  i 

u, l  max 

0.016.  This  mean  difference  Is  a  5.4  standard  deviation  result.  Using  this 
number,  we  can  compute  the  probability  that  the  18  Zeeman-like  frequency 
patterns  were  obtained  from  a  randomly  distributed  set  of  peaks  in 
frequency.  This  probability  is  less  than  10~^.  This  test  (which  is  a 
sufficient-condition  teat)  firther  strengthens  the  case  that  multiplets  have 
been  identified  in  the  power  spectra. 
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Fig  *  The  observed  rotational  splitting  linear  in  a  shown  as  function  of 
l  for  Intermediate  degree  f-modes  (Rabaey  and  Hill  1986).  Note 

that  the  rotational  splittings  are  negative.  The  point  plotted 
with  the  square  box  is  the  value  of  v'Qt  found  for  t  -  33  by  Hill 
(1986b). 


18 


(No.i/Nmox),'? 


Jt 


Fig  5  A  comparison  of  the  oBserved  ratio  N0t/Ninax  ^circle9^  and  the 
probability  p  (♦)  as  a  function  of  1  for  intermediate  degree  f- 
oodes  (Rabaey  and  Hill  1986).  The  quantity  NQjL  is  the  number  of 
peaks  classified  in  a  given  multlplet  and  N  is  the  maximum 
number  possible  for  that  multlplet.  The  probability  p  is  the 
chance  of  finding  one  or  more  peaks  from  a  random  distribution 
with  t  0.065  uHz  of  the  frequency  v0lB  given  by  Equation  5.  The 
solid  line  is  the  average  value  of  p  in  this  region  of  the 
frequency  spectrum. 
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1.4.4  Low-Degree  Gravity  Modes 

The  amplitudes  of  peaks  in  the  power  spectra  obtained  with  the 
differential  Doppler  studies  based  on  the  Crimean  observations  (Kotov  et  al. 
1983)  and  with  the  1979  differential  radius  observations  of  SCLERA  have  been 
examined  for  evidence  of  low-degree  gravity  modes  with  m  -  0 
eigenfrequencles  between  77.9  uHz  and  132.2  uHz.  In  this  analysis,  the 
families  of  peaks  in  the  respective  power  spectra  are  classified  as 
rotationally  split  oultiplets  when  (1)  they  exhibit  a  Zeeman-like  frequency 
pattern,  (2)  they  possess  the  minimum  quadratic  and  cubic  terms  in  m  in  the 
departure  from  uniform  spacing,  (3)  the  linear,  quadratic,  and  cubic  terms 
are  consistent  with  those  found  for  neighboring  multlplets  in  n  and  t,  (4) 
the  frequency  of  the  axlsymmetric  mode  is  consistent  with  the  frequencies  of 
corresponding  modes  in  neighboring  multlplets,  (5)  the  ratio  of  the 
differential  velocity  and  differential  radius  amplitudes  is  confined  to  a 
unique  value  for  each  value  of  l  (see  Section  1.3.5),  and  (6)  the 
classification  of  multlplets  identified  by  criteria  (l)-(5)  is  confirmed  by 
the  observed  horizontal  spatial  properties  being  characterized  by  m  (see 
Section  1.3.3). 

The  internal  rotation  of  the  Sun  implied  by  the  rotational  splitting 
results  of  Hill,  Bos,  and  Goode  (1982)  and  Hill  (1984b,  1985a)  was  used  as  a 
starting  point  in  the  outlined  search  for  identification  and  classification 
of  Internal  gravity  modes.  The  initial  value  assumed  for  the  rotational 
splitting  between  contiguous  modes  of  a  gravity  mode  multiplet  was  taken  as 
-2.9[1  -  1/1(1  *1)]  yHz.  This  approximate  expression  for  the  l-dependence 
is  based  on  the  theoretical  properties  found,  for  example,  in  the  work  of 
Gough  (1981).  Should  sets  of  peaks  be  found  with  a  self-consistent  set  of 
Zeeman-like  frequency  patterns  and  if  all  of  the  other  properties  considered 
in  the  search  with  regard  to  t  and  m  are  satisfied,  the  search  for  other 
solutions  need  not  be  carried  out.  This  is  because  the  observed  spatial 
properties  are  expected  to  establish  at  a  statistically  significant  level 
the  uniqueness  of  a  given  identification  and  classification. 


5.  This  is  a  summary  of  the  work  by  Hill  (1985b). 
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The  result  of  this  analysis  was  the  classification  of  31  gravity  node 
oultiplets  with  1  i  l  S  5.  Of  a  maximum  possible  number  of  235  modes 
belonging  to  the  31  multiplets,  152  were  classified.  These  multlplets  are 
Included  in  the  tabulations  of  Section  1.4.6. 

1.4.5  Intermediate-Degree  Gravity  Modes ^ 

The  first  evidence  of  the  excitation  of  intermediate-degree  solar 

gravity  modes  was  obtained  in  the  1970s  by  Hill  and  Caudell  (1979).  They 

found  evidence  which  indicated  the  excitation  of  a  number  of  the  gravity 

modes  with  frequencies  near  250  uHz  and  370  yHz  and  20  <  t  <  40.  These 

results  were  obtained  using  the  1973  solar  diameter  observations.  With  the 

availability  of  the  1979  differential  radius  observations,  with  their 

improved  signal- to- noise  ratio,  and  the  development  of  the  SCLERA  mode 

classification  program,  four  gravity  mode  multlplets  have  been  Identified 

with  v  .  -  350  uHz  and  of  intermediate  degree  (1  -  30).  The  preliminary 
n»  t 

results  have  been  reported  by  Hill  (1986c). 

The  criteria  used  in  the  identification  of  intermediate-degree  g-mode 
multlplets  were  quite  similar  to  those  used  by  Hill  (1984b)  in  the  study  of 
low-order,  low-degree  acoustic  modes.  The  principal  differences  have  to  do 
with  the  starting  values  of  v'  used  in  the  search  for  multiplets  and  the 

*•  1 1 

method  used  to  determine  the  m  -  0  mode  and  the  degree  of  the  multlplet. 

For  this  work,  the  values  of  v'  ,  obtained  by  Hill,  Bos,  and  Goode  (1982) 

n,i 

and  Hill  (1985b)  were  used  as  the  initial  values  assumed  in  a  search  for 
evidence  of  multiplets.  The  difference  with  regard  to  the  assignment  of  m 
and  t  arose  because  of  the  limited  number  of  multlplets  included  in  the 
current  analysis,  the  complexity  of  the  g-mode  spectrum,  and  the  reduced 
effectiveness  of  0^  (cf.  Section  1.3.3)  in  the  study  of  intermediate-degree 
g-modes.  Its  reduced  effectiveness  in  this  study  is  due  in  part  to  the 
closeness  of  v’  ,  to  4/d  (see  Hill,  Bos,  and  Goode  1982;  Hill  1985b)  and  in 
part  to  the  effect  of  atmospheric  seeing  on  the  for  m  >  20  (see  Section 
Ilia  of  Hill,  Alexander,  and  Caudell  1985). 


6.  Synopsis  of  work  by  Hill  (1986c) 
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The  technique  used  to  Identify  the  m  •  0  mode  and  to  determine  the 
value  of  l  for  a  given  multiplet  was  based  in  that  analysis  on  the 

properties  of  Y®,  where  Y®  is  the  average  of  Y®  over  the  detector  geometry 

used  by  Bos  (1982)  (see  Bos  and  Hill  1 983 ) •  as  discussed  in  Section  1.3*4. 

The  tentative  assignments  of  n,  t,  ra,  and  v  .  for  the  multi plets  for  which 

n  i  x* 

mode- coupling  evidence  was  found  are  given  in  Table  3-  The  standard  solar 
model  of  Saio  (1982)  was  used  in  the  assignments  of  n. 


Table  3 

Phase- locked  Gravity  Mode  Multlplets 
(tentative  assignments) 


Classification 

Vn,t 

(vHz) 

g14,28 

341.81 

g15,30 

341.43 

g13.30 

-  350 

0r  g15, 34 

(t  -  32) 

-  350 

1.4.6  Observationally  Derived  Eigenfrequency  Diagram 

It  is  quite  common  in  the  study  of  the  solar  5  min  oscillations  to 

display  the  eigenfrequencies  in  a  k  -  u  diagram.  In  practice,  this 

corresponds  to  plotting  the  eigenfrequency  as  a  function  of  the  wavenumber  k 

1  /2 

which  is  «  Ct(t*1)]  /R.  Such  information,  traditionally  based  on  Doppler 

shift  observations,  is  typically  confined  to  frequencies  between  2  mHz  and 
4  mHz. 

The  mode  classification  program  at  SCLERA  has  made  it  possible  to 
extend  the  observed  k  -  w  diagram  from  2  mHz  to  80  uHz  for  low-  and 
intermediate-degree  acoustic,  f- ,  and  gravity  modes.  This  mode 
classification  program  used  the  differential  radius  observations  from 
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SCLERA,  the  differential  velocity  and  infrared  observations  from  the  Crimean 
Astrophysical  Observatory  for  long  period  oscillations,  a  series  of  velocity 
observations  of  the  5  min  oscillations,  and  the  total  irradiance 
observations . 

The  results  of  this  mode  classification  program  are  summarized  in  Table 

4,  which  shows  the  mode  type,  order,  degree,  and  subsections  of  Section  1.4 

where  the  analysis  is  presented.  The  domain  of  the  results  are  also 

Indicated  in  Figure  6  In  a  \i.  *  l  diagram  analogous  to  the  k  -  u  diagram 

n ,  i 

discussed  above.  To  date,  166  multiplets  have  been  identified  and  •  1209 
resolved  modes  classified.  The  observed  m  -  0  eingenfrequencies  are  found 
to  be  -  10  pHz  above  the  theoretical  eigenfrequencies  for  the  acoustic  and 
f~modes  and  to  be  -  1  uHz  above  the  theoretical  eigenfrequencies  for  the 
low-degree  gravity  modes. 


Order 


Mode  Type 

Subsection 

n 

acoustic 

1.4.1 

1 2<n<27 

acoustic 

1.4.2 

1,2,3 

f 

1.4.3 

0 

g-avity 

1.4.4 

5<n<29 

«  •• 

gravity 

1.4.5 

n  »  1 5 

Table  4 


Degree 

l 

Number  of 
Multiplets 

Nunber  of 
Modes 

0<l<6 

83 

184 

2<l<22 

30 

293 

19<l<36 

18 

340 

1<l<5 

31 

152 

28<l<34 

4 

240 

1.5  Teats 

The  tests  of  the  SCLERA  observations  and  of  the  results  from  the 
multlplet  classification  program  fall  into  two  classes:  comparisons  of 
results  based  on  two  different  years  of  SCLERA  observations,  and  comparisons 
of  SCLERA  observations  with  those  based  on  observations  obtained  at 
different  observatories  using  different  observing  techniques. 

There  are  foir  sets  of  SCLERA  observations  that  have  been  analyzed 
during  the  last  five  years.  The  observations  are  those  made  in  1978  by 
Caudell  (1980),  in  1979  by  Bos  (1982),  in  1983  by  Ti  and  Czarnowskl  (1986), 
and  in  1985  by  Oglesby  (1986).  The  1979  observations  have  been  analyzed  by 
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The  theoretical  m  -  0  elgenfrequency  spectrum  based  on  the 
standard  solar  model  of  Saio  (1982).  The  enclosed  areas  represent 
the  domain  of  tne  mode  classifications  that  have  been  obtained  in 
the  SCLERA  mode  classification  program  (see  Table  H) .  The 
acoustic  and  yavity  modes  are  indicated  by  positive  and  negative 
radial  order  numbers,  respectively. 


Hill  (1984b)  and  30  low-order,  low-degree  acoustic  mode  oultiplets 
classified  (cf.  Section  1.4.2).  The  1978,  1983,  and  1985  observations  have 
been  compared  to  the  1979  observations  by  testing  for  signals  in  these 
former  three  sets  of  observations  at  the  eigenfrequencies  derived  from 

the  1979  observations .  Statistically  significant  positive  results  have  been 
obtained  in  each  of  the  three  comparisons.  These  findings,  reviewed  in 
Sections  1.5.1.  1.5.2,  and  1.5.3,  respectively,  are  the  results  of 
particularly  stringent  tests  of  the  quality  of  the  data  and  the  accuracy  of 
the  mode  classification  program. 

Results  based  on  the  1979  SCLERA  observations  are  also  compared  with 
results  obtained  elsewhere:  with  infrared  observations  made  at  the  Crimean 
Astrophysical  Observatory,  with  the  total  irradiance  observations  made  on 
the  Solar  Maximum  Mission  (SMM) ,  and  with  a  series  of  velocity  observations 
of  the  5  min  oscillations.  Each  of  these  comparisons,  also  reviewed  below, 
leads  to  positive  results. 

7 

1.5.1  Comparison  of  1978  and  1979  SCLERA  Observations 


The  1978  solar  diameter  observations  of  Caudell  et  al.  (1980)  have  been 
examined  for  evidence  of  the  low-order,  low-degree  acoustic  modes  found  by 
Hill  (1984b)  in  the  1979  observations  of  Bos  and  Hill  (1983).  The  power 
spectrum  of  the  1978  observations  was  tested  for  evidence  of  (1)  a  non- 
random  distribution  of  peaks  relative  to  the  eigenfrequency  spectrum  found 
by  Hill  (1984b),  (2)  spatial  symmetry  properties  of  the  eigenfunctions  for 
the  modes  classified  by  Hill  (1984b),  and  (3)  a  e-dependence,  specified  by 
angular  order  m,  of  the  amplitudes  of  oscillations,  a,  for  these  same  modes. 
The  character! zation  of  the  symmetry  properties  of  the  eigenfunctions  and 
the  values  of  m  were  taken  from  the  work  of  Hill  (1984b).  The  findings  of 
these  three  tests  confirm  the  detection  of  modes  and  their  classifications 
as  given  by  Hill  (1984b). 

The  first  step  in  the  comparison  of  the  1978  and  1979  observations  was 
to  ascertain  whether  the  power  spectrum  for  the  independent  set  of  1978 


7.  This  is  a  synopels  of  the  work  by  Hill  and  Caudell  (1985) 
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observations  consisted  of  a  set  of  peaks  randomly  distributed  In  frequency 

relative  to  the  multiplet  eigenfrequencies  found  by  Hill  (1984b).  The 

negative  result  obtained  represents  a  confirmation  at  the  2.7  o  level  of  the 

multiplet  classifications  of  Hill  (  1984b)  (cf.  Sections  IIIc  and  Hid  of 

Hill  and  Caudell  1985).  The  second  step  was  to  determine  whether  the 

spectrum  for  the  1978  observations  was  consistent  with  the  symmetry 

properties  expected  for  the  eigenfunctions  of  the  multiplet  modes  classified 

by  Hill  (1984b).  The  spectrum  was  found  to  be  consistent  at  the  4.1  <j 

level.  A  third  independent  test  substantiates  the  isomorphism  found  by  Hill 

(1984b)  between  the  observed  square  amplitude  of  the  modes  and  |m|.  The 

results  of  this  test  support  the  above  confirmations.  Further,  the  positive 

2 

cross-correlation  coefficient  obtained  between  the  respective  <a  > 
establishes,  at  a  1.8  a  level,  further  support  for  the  conclusions  by  Hill 
(1984b)  regarding  multiplet  detection  and  classification. 

These  three  results,  one  at  a  confidence  level  of  2.7  o,  a  second  at 
4.1  a,  and  a  third  at  1.8  o,  are  all  independent.  This  is  apparent  from  the 
fact  that,  in  one  case,  a  search  was  made  for  peaks  at  frequency  locations 
where  they  can  be  present  on  symmetry  grounds  and  in  the  second  case,  a 
search  was  made  for  peaks  at  frequency  locations  where  they  cannot  be 
present  on  symmetry  grounds.  In  the  third  case,  a  search  was  made  for  an 
isomorphism  between  the  detected  amplitude  and  |m|.  Thus,  when  combined, 
these  three  independent  results  present  a  very  strong  test  of  the  internal 
consistency  of  the  multiplet  classification  program  with  positive  results. 

O 

1.5.2  Comparison  of  1983  and  1979  SCLERA  Observations 

The  1983  data  was  obtained  with  the  same  detector  used  in  1979  at 
SCLERA.  However,  in  1983  the  detector  was  rotated  periodically  so  that  the 
edge  of  the  solar  disk  could  be  sampled  at  eight  different  positions.  These 
eight  positions  were  chosen  to  lie  on  four  solar  diameters  which  were 
oriented  at  -90,  -45,  0,  and  45  degrees  from  the  projected  solar  axis. 


8.  This  is  a  synopsis  of  the  work  by  Yi  (1986).  Preliminary  results  have 
been  reported  in  Yi  and  Czarnowski  (1986). 
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The  1983  data  spans  52  days,  from  30  May  1983  to  20  July  1983.  The 
actual  analysis  used  a  total  of  319.5  hours  of  data  drawn  from  44  days. 

The  detector  was  rotated  to  a  new  diameter  every  64  seconds.  At  each 
diameter,  the  solar  edge  was  scanned  8  times.  Since  the  detector  has  6 
slits  arranged  in  a  pattern  that  is  symmetric  with  respect  to  the  solar 
diameter  (see  Bos  and  Hill  1983),  48  limbs  were  sampled  and  recorded  each 
time  the  detector  was  rotated  to  a  new  position.  Thus,  for  the  equatorial 
position  (and  each  of  the  other  positions)  8  limbs  were  recorded  every  256 
seconds. 

The  techniques  used  at  SCLERA  to  obtain  the  FFTD  lock-on  points  for 
each  limb  have  been  extensively  discussed  elsewhere  (cf.  Bos  and  Hill  1983). 
The  spectra  analyzed  were  obtained  with  a  standard  FFT  program  applied  to 
the  lock-on  points. 

Previous  work  at  SCLERA  (cf.  Hill  1984b)  produced  a  set  of  solar 
acoustic  mode  eigenfrequencies  based  on  the  analysis  of  the  1979  data.  The 
1983  data  were  examined  for  evidence  of  these  modes  in  order  to  further  test 
the  mode  classification  program  and,  should  a  positive  result  be  obtained, 
possibly  learn  something  about  stability  of  the  frequencies  of  the  low- 
order,  low-degree  acoustic  modes  over  a  five-year  period. 

In  this  analysis,  the  previously  blassified  n  *  1  rauTtiplets  with 

l  -  4 . 22  were  divided  into  three  groups:  1)  l  •  4 . 9;  2)  l  - 

1 . . 16;  and  3)  i  -  17,...,  22.  It  has  been  found  that  slight  shifts  of 

the  eigenfrequencies  based  on  the  1979  observations  match  at  a  statistically 
significant  level  a  set  of  peaks  in  the  1983  spectra.  The  shifts  in 
frequencies  which  were  found  effective  are  -0.02,  0.03,  and  0.05  uHz  for 
groups  1,  2,  and  3,  respectively.  A  tolerance  of  ±0.06  uHz  was  used  to 
determine  whether  a  peak  in  a  spectrum  matched  a  given  eigenfrequency. 

A  second  comparison  of  the  1983  and  1979  observations  was  made  based  on 
an  assumed  (exp  im^)- dependence  of  the  eigenfunctions.  The  information  on 
the  ^-dependence  is  contained  in  D^j  ,  which  is  a  function  of  m  (cf.  Section 
1.3.3;  Hill  1984b;  Hill,  Alexander,  and  Caudell  1985).  These  were  divided 
into  three  groups  where  |m|  -  0,...,7;  |m|  -  8,.. .,14;  and  |m|  -  15,. ..,22 
and  the  weighted  average  <Dij>  ar®  0.125  ±  0.043,  -0.0079  ±  0.074,  and  0.476 
±  0.034  for  the  three  groups,  respectively.  The  <D^>'s  obtained  for  the 
three  groups  are  consistent  at  a  statistically  significant  level  with  the 
corresponding  quantities  calculated  using  the  1979  observations. 
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Also,  In  previous  work  at  SCLERA  (cf.  Hill  1984b),  It  has  been  shown 
that  another  quantity  —  R  —  can  be  used  In  determining  the  orientation 

•  J 

of  an  axis  of  symmetry  in  the  Sun.  Use  was  made  of  this  quantity  in 

analyzing  the  1983  data.  A  total  of  182  acoustic  modes  with  1  -  4 . 22 

were  Identified  in  the  1983  data.  Those  with  1-8,  9,  10,  11,  12  were  used 
to  calculate  the  orientation  of  the  internal  axis  of  symmetry.  A  total  of 
29  modes  with  these  1  values  were  found.  The  angle  of  the  axis  of  symmetry 
from  the  pole,  calculated  from  the  Rjj's  observed  for  these  modes,  is 
0.2  t  0.8  deg. 

1.5.3  Comparison  of  1985  and  1979  SCLERA  Observations^ 

SCLERA  has  been  Involved  since  1984  in  expanding  its  differential 
radius  observational  program  with  the  introduction  of  a  new  technique  to 
search  for  global  solar  oscillations  in  the  visible  to  near  Infrared 
continuum  (cf.  Section  1.7).  In  1985,  data  was  collected  using  this  new 
technique. 

Two  hundred  ninety-one  hours  of  data  were  collected  in  a  70-day  period 
from  23  March  to  31  May  1985.  A  200  x  200  arcsec  aperture  was  centered  on 
the  solar  disk  and  the  Intensity  at  48  different  wavelengths  was  measured 
within  the  spectral  range  of  0.5  to  1.7  u.  The  radiation  intensities 
observed  at  these  wavelengths  were  Fourier  analyzed  and  the  resulting  power 
spectra  examined  for  evidence  of  the  low-order  ,  low-degree  acoustic  modes 
identified  in  the  1979  SCLERA  observations  by  Hill  ( 1 984b ) . 

The  eigenfrequency  spectrum  for  n-1,  l  -7,. ..,22  multi plets  was 
divided  into  the  three  groups  7S  l  S  12,  13  S  i  S  17,  and  18  i  l  i  22.  The 
elgenfrequencies  of  the  modes  belonging  to  a  group  were  shifted  in  frequency 
as  a  group  to  search  for  a  maximum  coincidence  rate  between  the  peaks  in  the 
power  spectra  of  the  1985  data  and  the  frequencies  of  the  previously 
classified  modes.  A  maximum  allowable  shift  was  ±0.10  uHz.  Maxima  in  the 
coincidence  rates  were  fourd  for  the  frequency  shifts  of  -0.06,  0.01,  and 
0.04  uHz,  respectively. 


9.  This  is  a  synopsis  of  work  by  Oglesby  (1986).  Preliminary  results  have 
been  reported  in  Oglesby  (1986). 


27 


The  nean  difference  d  between  the  observed  coincidence  rate  and  the 
expected  rate  between  two  distributions  that  are  randomly  distributed  with 
respect  to  each  other  was  calculated.  This  yielded  a  value  for  c  for  each 
of  the  three  groups  of  multi plets  of 

0.214  t  0.046;  7  i  t  i  12 

d  -  0.145  ±  0.040;  13  S  l  i  17  .  (6) 

0.132  t  0.035;  18  i  t  *  22 

These  values  of  d  are  4.7,  3.6,  and  3.8  standard  deviations,  respectively, 
above  zero.  The  probability  of  obtaining  a  particular  positive  3/ o  where 

the  v  .  of  the  classified  modes  for  the  three  bins  has  been  allowed  to  vary 

ni 

by  ±0.10  uHz  is  estimated  to  be  within  an  order  of  magnitude  of 

1  x  10"5;  7  S  l  S  12 

p  <  2  x  1<f3  ;  13  i  l  S  17  .  (7) 

4  x  10~4  ;  18  $  l  *  22 

For  the  three  groups  of  multlplets  described  above,  the  positive 
results  at  the  4.7o,  3.6o,  and  3>8o  level  represent  a  highly  statistically 
significant  result  which  leads  to  a  rejection  of  the  hypothesis  that  the 
peaks  in  the  power  spectrum  of  the  1985  data  are  randomly  distributed  with 
respect  to  the  classified  eigenfrequency  spectrum  based  on  the  1979 
observations.  This  represents  a  particularly  important  finding  since  the 
confirmation  comes  from  an  independent  set  of  observations  using  a  different 
observational  technique. 

It  was  also  noted  in  the  analysis  that  apparently  the  eigenfrequencies 
underwent  a  slight  frequency  shift  dependent  on  1  between  1979  and  1985. 
Since  the  multlplets  discussed  here  reside  primarily  in  the  convection  zone, 
knowledge  of  any  change  of  these  eigenfrequencies  may  yield  interesting 
information  about  the  level  of  variations  in  the  properties  of  the 
convection  zone  over,  for  example,  the  solar  cycle. 

Recently,  f-modes  have  seen  identified  in  the  1979  data  (cf.  Section 
1.4.3).  An  analysis  similar  to  the  one  above  is  underway  to  determine  if  a 
correlation  for  these  mooes  exists  between  1979  and  1985.  Preliminary 
results  Indicate  agreement  at  the  3o  level. 
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1.5.**  Comparison  of  Infrared  and  Differential  Radius  Observations10 

The  well-documented  160  min  period  solar  oscillation  (Kotov  et  al . 
1983)  has  been  classified  by  Hill.  Tash,  and  Padln  (1986)  as  an  internal 
gravity  mode  with  l  -  2  and  m  -  2  where  t  is  the  degree  and  m  the  angular 
order  of  the  eigenfunction.  This  particular  classification  was  determined 
by  comparing  the  results  obtained  on  the  160  min  period  oscillation  from 
differential  velocity  observations  made  at  the  Crimean  Astrophys  1  cal 
Observatory  and  from  differential  radius  observations  made  at  SCLERA.  The  t 
and  m  classification  has  subsequently  been  confirmed  in  the  work  of  Hill 
(1985b),  in^  which  the  160  min  period  oscillation  is  further  classified  as 
the  m  -  2  member  (retrograde)  of  the  g  multlplet. 

Differential  infrared  observations  of  the  160  min  period  oscillation 
were  made  at  the  Crimean  Astrophysical  Observatory  in  1977  and  1978  by  Kotov 
and  Koutchmy  (1979)  and  in  1981  by  Kotov,  Koutchoy,  and  Koutchmy  (1983). 
The  above  m  classification  for  this  mode  was  tested  further  by  comparing 
these  Infrared  observations  with  the  differential  '•adius  observations  made 
at  SCLERA  in  1979.  Such  an  analysis  Involves  comparison  of  relative 
amplitudes  and  phases  of  signals  obtained  at  the  two  observatories .  Besides 
furnishing  an  additional  test  to  the  g  classification,  a  comparison 
like  this  presents  the  opportunity  to  study  further  the  properties  of  the 
160  min  period  temperature  eigenfunction  in  the  photosphere. 

For  this  comparison,  predicted  Infrared  slpials  based  on  the  SCLERA 
1 979  observations  are  compared  to  those  observed,  assuming  |m|  -  2,  the 
value  assigned  to  the  160  min  period  oscillation  by  Hill,  Tash,  and  Padln 
(1986).  The  amplitude  and  phase  of  the  predicted  160  mir.  period  signal  for 
the  1977-1978  observations  are  in  good  agreement  with  the  signal  observed. 
The  predicted  amplitude  agrees  with  the  amplitude  found  in  the  1981 
observations;  however,  no  such  agreement  is  found  for  the  relative  phases. 
On  the  other  hand,  the  1977-1978  infrared  observations  are  of  a  much  higher 
statistical  quality  than  are  the  1981  observations.  If  the  1977-1978  and 
1981  observations  were  weighted  according  to  observing  time  alone  (leaving 
out  the  Important  difference  in  frequency  resolution),  the  relative  weights 


10.  This  is  a  synopsis  of  the  work  by  Hill  (1986  a). 
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would  be  4.9  to  1.  Taking  this  weighting  into  account,  it  is  concluded  that 
both  amplitudes  and  phases  of  observed  infrared  signals  are  in  good 
ag-eeaent  with  those  predicted.  Furthermore,  the  relative  phases  of  the 
1977-78  Infrared  observations  and  the  1979  differential  radius  observations 
are  consistent  with  respect  to  phase  only  if  |s|  •  2.  Vie  thus  have  an 
Independent  confirmation  of  the  |m|  •  2  assignment  for  the  160  min  period 
oscillation  made  by  Hill,  Tash,  and  Padin  (1986). 

1.5.5  Comparison  of  Total  Irradiance  and  SCLERA  Observations11 

A  power  spectrum  of  total  Irradiance  observations  has  been  presented  by 
Woodard  and  Hudson  (1983)  in  the  frequency  range  of  gravity  modes.  For  the 
same  frequency  region,  a  number  of  gravity  modes  have  been  classified  by 
Hill  (1985b)  (cf.  Section  1.4. A)  using  differential  radius  observations  and 
the  differential  velocity  observations  of  Kotov  et  al.  (1983).  A  comparison 
has  been  made  between  the  total  irradiance  power  spectrum  and  the  classified 
gravity  mode  spectrum  to  test  the  mode  classification  program  and  to  test 
for  evidence  of  gravity  modes  in  the  total  irradiance  observations.  In  a 
second  test  for  such  evidence,  a  comparison  is  made  between  the  amplitudes 
of  signals  in  the  total  irradiance  and  differential  radius  observations. 
The  first  test  represents  a  study  of  the  frequency  properties  of  the 
observed  sigials  while  the  second  test  represents,  because  of  differences  in 
the  spatial  filter  fmetions  of  the  two  different  types  of  observations ,  a 
study  of  the  horizontal  spatial  properties  of  the  observed  signals. 

The  comparison  of  the  classified  g*avlty  mode  spectrum  with  the  total 
irradiance  power  spectrum  entailed  the  examination  or  the  peaks  in  the  total 
irradiance  power  spectrum  and  the  frequencies  of  the  classified  gravity 
modes.  It  was  observed  that  no  two  or  more  members  of  the  set  of  gravity 
mode  frequencies  have  a  common  nearest  peak  in  the  total  Irradiance 
spectrum,  and  in  10  of  the  11  cases,  the  magnitude  of  the  difference  between 
the  gravity  mode  frequencies  and  the  total  irradiance  peak  frequencies  was 
<  0.08  uHz.  The  root-man- square  deviation  was  0.065  uHz. 


11.  This  is  a  synopsis  of  the  work  by  Hill  and  Kroll  (1986). 
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The  statistical  significance  of  the  frequency  association  described 
above  was  determined  by  estimating  the  probability  of  obtaining  the  saae 
association  between  two  unrelated  sets  of  computer-generated  frequencies. 
Based  on  162, 055  sets  of  such  frequencies,  it  is  inferred  that  there  is  a 
probability  of 


p  -  0.020  (8) 

for  obtaining  the  observed  frequency  alignment  and  association  between  the 
total  irradlance  spectrum  and  the  relevant  portion  of  Hill's  (1985b) 
classified  gravity  mode  spectrum  if  the  two  spectra  are  unrelated.  This 
small  percentage  is  evidence  in  support  of  the  correctness  of  the  gravity 
mode  classification  of  Hill  (1985b)  and  of  the  importance  of  the  effects  of 
gravity  modes  in  the  total  irradlance  spectrum. 

In  general,  the  ratio  of  amplitudes  of  an  oscillation  obtained  in  two 
different  types  of  observations  will  depend  on  both  1  and  m  of  the  mode, 
principally  because  of  different  detector  spatial  filter  functions.  Thus, 
if  the  relevant  modes  of  oscillations  have  not  been  classified,  the 
comparison  of  two  different  types  of  observations  can  be  very  difficult. 
However,  if  the  g*avlty  mode  spectrum  classified  by  Hill  (1985b)  is  correct, 
it  can  be  used  to  reduce  the  ratio  of  amplitudes  obtained  in  two  different 
types  of  observations  from  a  two- parameter  function  to  a  one-parameter 
function.  9uch  a  reduction  in  independent  variables  could  make  it  possible 
to  test  further  for  the  correctness  of  the  mode  classification  prop*am  and 
for  the  presence  of  gravity  mode  signals  in  the  total  irradlance 
observations. 

The  ratios  of  the  amplitudes  of  signals  for  total  irradlance 
observations  and  differential  radius  observations  were  examined  for 
agreement  with  those  expected  based  on  the  mode  classifications  of  Hill 
(1985b).  It  was  found  l n  this  analysis  that  the  probability  of  obtaining 
the  observed  correlation  of  -elative  amplitudes  is 

P  v  «  10' 3  (9) 

if  the  mode  classif icatio-s  -ere  incorrect  and/or  if  the  total  irradlance 
power  spectrum  contains  no  evidence  of  p*avity  modes. 
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The  results  of  the  second  test,  combined  with  the  Independent  findings 
of  the  first  test,  furnish  a  quite  statistically  significant  result  in  terms 
of  the  correctness  of  the  gravity  mode  classification  In  both  l  and  m  by 
Hill  (1985b)  and  a  quite  statistically  significant  demonstration  that 
evidence  of  gravity  modes  is  present  in  the  total  irradlance  power  spectrum. 

12 

1.5.6  Test  of  e  Uniformly  Rotating  Convection  Zone  Model. 

Considerable  discrepancy  exists  among  the  reported  determinations  of 
rotational  splitting  in  the  eigenfrequency  spectrum  of  the  Sun.  The  values 
obtained  by  Claverle  et  al.  (1981)  and  Woodard  and  Hudson  (1983),  and  their 
inconsistency  with  the  results  of  Hill,  Bos,  and  Goode  (1982)  and  Hill 
(1985a),  are  discussed  in  Hill  (1985a).  The  sharpest  discrepancy  is  found 
between  the  results  of  Duvall  and  Harvey  (1984),  Brown  (1985),  and  Llbbrecht 
(1986),  obtained  for  5  min  oscillations,  and  those  of  Hill  (1984,  1985a, b), 
obtained  for  5  min,  •  27  min,  and  •  160  min  oscillations.  The  works  of 
Duvall  and  Harvey  (1984),  Brown  (1985),  and  Llbbrecht  (1986)  use  Doppler 
shift  observations,  while  the  three  works  of  Hill  use  differential  velocity 
observations,  made  at  the  Crimean  Astrophysical  Observatory  (Kotov  et  al. 
1  983).  and  differential  radius  observations  made  at  SCLERA.  The  resolution 
of  this  disagreement  must  depend  upon  an  understanding  of  the  complex 
details  of  the  respective  data  acquisition  systems  and  the  assumptions  and 
procedures  used  in  analysis.  Although  this  discrepancy  has  not  yet  been 
resolved,  a  test  has  been  made  by  Hill,  Roeenwald,  and  Rabaey  (1985)  and 
extended  by  Hill  (1986b)  of  the  model  that  the  internal  rotation  of  the 
convection  zone  is  correctly  given  by  the  rotation  curve  obtained  by  Duvall 
et  al.  (1984). 

It  is  apparent  from  the  properties  of  the  splitting  kernels  that  the 
Internal  rotation  curve  obtained  by  Duvall  et  al.  (1984)  predicts  a  synodic 
rotational  splitting  of  -  0.425  uHz  for  the  low-order,  low-deg-ee  acoustic 
mode  and  intermedlate-deg*ee  f-mode  multiplets  with  eigenfrequencies  between 
450  uHz  and  600  uHz.  This  featvre  of  the  Duvall  et  al.  (1984)  curve,  which 
is  oommon  to  all  essentially  uniformly  rotating  models  of  the  convection 


12.  This  is  a  synopsis  of  the  work  by  Hill  (1986b). 
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zone,  makes  it  particularly  easy  to  test  for  the  correctness  of  this  class 
of  models.  Note  that,  for  a  sizable  fraction  of  the  observed  signals,  the 
frequency  spacings  between  power  spectra  peaks  due  to  these  oscillations 
should  be  either  mod  -  0.  425  yHz  or  mod  -  0.85  uHz,  depending  upon  tne 
spatial  filtering  properties  of  the  signal  detector.  This  follows  as  a 
direct  consequence  of  the  fact  that,  for  oscillations  whose  coherence  times 
are  longer  than  the  observing  time  span,  the  shape  of  the  power  spectrum  for 
each  oscillation  detected  is  determined  completely  by  the  observing  window 
function,  i.e.,  it  is  the  same  for  each  of  these  individual  modes  of 
oscillation.  A  mod  -  0.425  uHz  or  mod  -  0.85  uHz  characteristic  peak 
spacing  is  a  relatively  easy  featire  to  identify. 

The  results  of  the  search  for  evidence  of  mi  form  ly  rotating  convection 
zone  models  are  negative.  At  the  10.6  a  level,  no  evidence  for  an 
essentially  uniformly  rotating  convection  zone  was  found  if 

<Q  (r)/2*>  i  0.425  uHz  ,  (10) 

eq 

and  at  the  4.8  o  level,  no  evidence  was  found  for  an  essentially  uniformly 
rotating  convection  zone  if 

<Q  (r  )/2w>  S  0.438  uHz  (11) 

eq 

where  is  the  rotational  speed  at  the  equator. 

1.6  Comparison  of  Results  from  Different  Types  of  Observations: 

Work  on  Resolution  of  Problems 

One  of  the  more  Important  activities  at  SCLERA  has  been  the  search  to 
understand  some  of  the  problems  that  have  arisen  in  comparisons  of  results 
from  different  observations.  The  problems  addressed  to  date  concern  the 
differences  in  reported  multiplet  fine  structure  and  in  inferred  internal 
rotation,  the  varied  results  on  the  160  min  oscillation,  the  relatively 
large  scatter  in  the  reported  elgenfrequencies  of  the  low-deg-ee  five-mln 
oscillations,  and  the  reported  detection  of  intermediate-degree  gravity 
nodes.  This  work,  which  is  siamarized  in  the  following  sections,  has  served 
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to  more  clearly  define  the  features  of  several  of  the  problems  and  to 
introduce  interesting  hypotheses  for  some  others. 

1.6.1  Multlpleb  Fine  Structure  and  Inferred  Internal  Rotation 

The  Identification  of  the  multi  pi et  fine  structure  due  to  slow  internal 
rotation  of  the  Sun  has  presented  many  more  problems  than  anticipated,  and 
as  a  result  the  analysis  technique  of  Section  1.3.1  is  not  as  easy  to 
implement  as  once  thought.  In  particular,  there  is  widespread  observational 
disagreement  among  the  reported  determinations  of  rotational  splitting  in 
which  there  are  essentially  three  different  sets  of  findings  (see  Hill,  Bos, 
and  Goode  1982,  Hill  1984b,  1985a, b,  and  Rabaey  and  Hill  1986  for  one 
finding;  Duvall  and  Harvey  198**,  Brown  1985,  and  Libbrecht  and  Zirin  1986 
for  a  second;  and  Duvall,  Harvey,  and  Pomeranta  1986  for  a  third).  The 
first  set  of  findings  is  essentially  based  on  the  differential  radius 
observations  of  typically  longer  period  oscillations,  the  second  set  of 
findings  is  based  on  Doppler  shift  observations  of  the  five  min 
oscillations,  and  the  third  set  of  findings  is  based  on  intensity 
observations  within  an  absorption  line  of  the  five  min  oscillations. 

The  dissonance  of  the  Information  about  rotational  'effects  and  the 
apparent  Rouping  of  the  different  findings  according  to  observing  technique 
demonstrates  the  merits  of  employing  different  observational  techniques  in 
testing  the  validity  of  a  given  inter pretation  of  a  set  of  observations. 
The  need  for  continued  diversity  in  observational  technique  for  future 
studies  is  vital  for  the  resolution  of  this  difficulty  and  for  the  further 
development  of  solar  seismology.  It  is  in  this  spirit  that  a  number  of 
different  analyses  have  been  undertaken  with  regard  to  the  Identification  of 
the  proper  rotational  splitting  effects. 

The  analyses  fall  into  several  categories.  They  are:  (1)  the 
determination  of  the  multlplet  fine  structure  for  low-des"ee  five  min  modes, 
low-degree  acoustic  modes,  intermediate-degree  f-modes,  and  low-  and 
intermediate- degree  gravity  modes  using  primarily  the  differential  radius 
observations  from  SCLERA  based  on  the  SCLERA  mode  classification  program; 
(2)  the  incorporation  of  the  observations  obtained  using  different  observing 
techniques  into  the  SCLERA  mode  classification  program;  (3)  the  inferring 
from  these  results  the  implied  internal  rotation  and  the  determination  of 
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whether  the  observed  multlplet  fine  structure  is  consistent  with  a  single 
internal  rotation  curve  for  the  Sun,  (4)  the  comparison  of  the  results 
obtained  from  one  year's  observation  with  those  of  another  as  a  check  on  the 
mode  classification  program,  and  (5)  the  introduction  of  a  new  observational 
technique  to  enlarge  the  number  of  different  types  of  observations  used  to 
study  this  problem. 

The  results  from  the  works  that  fall  into  these  five  categories  have 
not  led  to  a  resolution  of  the  rotational  splitting  problem  but  have  served 
to  more  clearly  define  the  nature  and  extent  of  the  difficulty. 

The  multiplet  fine  structure  for  the  low-order,  low-degree  acoustic 
modes,  intermediate-degree  f-modes,  and  intermediate-degree  gravity  modes 
obtained  in  the  analysis  of  the  1979  differential  radius  observations  is 
taken  from  Hill  (1984b),  Rabaey  and  Hill  (1986),  and  Hill  (1986c).  The 
m  -  0  eigenfrequency  for  the  modes  is  discussed  in  Sections  1.4.2,  1.4.3, 
and  1.4.5,  respectively. 

The  multiplet  fine  structure  for  the  five  min  modes  has  been  obtained 
in  an  analysis  based  on  differential  radius  observations  combined  with  a 
number  of  Doppler  shift  observations  and  total  irradiance  observations 
(Woodard  and  Hudson  1983,  Claverie  et  al.  1981,  Grec,  Fossat,  and  Pomerantz 
1983,  Duvall  and  Harvey  1983,  and  Scherrer  et  al.  1983). 

The  multiplet  fine  structure  for  the  low-degree  gravity  modes  has  been 
obtained  in  an  analysis  based  on  differential  radius  observations  combined 
with  the  differential  velocity  observations  made  at  the  Crimean 
Astrophyslcal  Observatory  (Kotov  et  al.  1983).  The  basis  of  this  analysis 
is  discussed  in  Section  1.6.2. 

The  degree  of  internal  consistency  of  the  multiplet  fine  structure  has 
been  examined  in  a  series  of  works  (Hill  et  al.  1984,  Hill,  Rosenwald,  and 
Rabaey  1985,  and  Hill  et^  a_l.  1986),  in  which  each  succeeding  analysis 
included  the  input  of  a  larger  number  of  multi plets.  In  the  most  recent 
work  (Hill  et  al.  1986),  the  multiplet  fine  structure  for  30  low-order,  low- 
degree  acoustic  modes,  for  83  low-degree  5  min  multiplets,  and  for  two 
intermediate-degree  gravity  mode  multiplets  was  used.  These  analyses 
indicated  that  the  multiplet  fine  structure  obtained  for  the  multiplets  is 
consistent  with  a  single  internal  rotation  curve  for  the  Sun.  They  also 
indicate  that  the  multiplet  fine  structure  obtained  for  the  intermediate- 
degree  f-modes  is  consistent  with  this  internal  rotation  curve. 
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The  efficiency  of  the  multiplet  classification  program  that  yielded  the 
multiplet  fine  structure  has  also  been  tested  by  comparing  the  power  spectra 
of  observations  obtained  in  different  years  and  with  different  observing 
techniques.  These  tests  are  described  in  Sections  1.5.1,  1.5.2,  1.5.3,  and 
1.5.5.  It  is  the  first  three  of  these  four  tests  that  furnish  the  most 
critical  test  of  the  multiplet  classifications.  This  is  because  these  tests 
are  concerned  with  the  low-order  ,  low-degree  acoustic  modes  which  have  a 
coherence  time  <  1  y.  Consequently,  these  tests  are  not  seeking  to 
determine  if  one  power  spectrum  is  random  or  not  random  relative  to  another 
power  spectrum  (which  would  not  test  the  mode  classifications),  but  are 
trying  to  ascertain  if  a  power  spectrum  is  random  or  not  random  relative  to 
the  frequency  spectrum  of  the  classified  modes  based  on  the  1979 
observations.  As  noted  in  Sections  1.5.1,  1.5.2,  and  1.5.3,  positive 
results  were  obtained  in  each  of  the  three  tests,  indicating  a  high  level  of 
proficiency  for  the  mode  classification  program. 

These  series  of  analyses  and  tests  have  yielded  an  observationally 
based  multiplet  fine  structure  that  is  internally  self-consistent.  Further, 
the  structure  is  observed  to  be  stable  from  one  year  to  another.  There  is 
one  common  feature  of  the  observations  used  in  this  analysis  that  may  be 
relevant  to  defining  the  nature  of  the  reported  differences  in  the  multiplet 
fine  structure.  That  feature  is  that  the  differential  radius  observations, 
the  differential  velocity  observations  (see  Sections  1.4.4  and  1.6.2),  and 
the  total  irradiance  observations  are  sensitive  primarily  to  the  temperature 
eigenfunctions  of  the  oscillations. 

1.6.2  Apparent-Velocity  Hypothesis!^ 

Comparing  observations  of  oscillations  detected  by  Doppler  shift 
measurements  with  ones  detected  by  solar  diameter  measurements  was  first 
thought  to  Involve  only  S„/r,  the  fractional  change  in  the  radius  and  its 
accompanying  horizontal  displacement,  ,  produced  by  an  oscillation.  The 
only  complexity  in  this  case  would  arise  from  the  radial  dependence  of  ^,/r 


13.  This  is  taken  from  the  wo r*  of  Hill,  Tash  and  Padln  (1985). 


and  £  /r  and  the  angular  dependence  of  the  oscillation  given  by  the 
n 

spherical  harmonic 

It  is  now  clear,  however,  that  another  dimension  introduced  further 
complexities.  Examination  of  the  observed  amplitudes,  from  both  Doppler 
shift  and  differential  radius  measurements,  reveals  yet  another  significant 
internal  discrepancy.  As  noted  above,  the  first  such  discrepancy  led  to  the 
recognition  that  the  solar  diameter  measurements  were  detecting  only  shape 
changes  in  the  limb- darkening  function.  The  second  major  discrepancy  that 
has  been  identified  by  Hill  and  Tash  (1984)  and  Hill,  Tash,  and  Padin  (1986) 
between  observed  amplitudes  and  theory  led  to  a  reexamination  of  the 
physical  mechanism  by  which  the  Doppler  shift  measurements  detect  the  long- 
period  oscillations.  It  was  quickly  discovered  in  this  reexamination  that 
the  apparent  Doppler  shift,  due  to  the  contribution  of  a  perturbation  in  the 
radiation  intensity  found  in  the  differential  radius  measurements  and  the 
equatorial  surface  velocity,  was  large  enough  to  give  rise  to  apparent 
velocities  of  fractions  of  m/sec  for  the  nonradlal  oscillations,  velocities 
of  the  same  order  as  those  reported  in  Doppler  shift  studies  (Hill,  Tash, 
and  Padin  1986). 

The  hypothesis  has  been  examined  that  the  primary  physical  mechanism  by 
which  the  1 60  min  long-period  solar  oscillation  is  detected  via  Doppler 
shifts  Involves  temperature  pertur batlons  conjoined  with  solar  surface 
rotation,  and  not  physical  displacement  of  the  solar  surface.  The  formalism 
has  been  developed  by  Hill,  Tash  and  Padin  (1986)  for  computing  the  apparent 
Doppler  shifts  from  the  temperature  pertir  bations  obtained  in  diameter  and 
differential  radius  observations.  Calculations  take  into  account  opacity 
effects  based  on  a  non^LTE  analysis  of  the  Fel  5124  transition  used  in  the 
differential  velocity  observations.  A  symmetric  line  profile  is  not 
assumed.  The  predicted  apparent  velocities  are  obtained  with  uncertainties 
of  +  120%,  -40%,  and,  for  physically  plausible  values  of  i,  the  predicted 
values  are  consistent  with  observed  velocities.  The  observed  velocities  for 
the  Crimean  and  Stanford  observations  of  the  1 60  min  oscillation  are  0.58 
m/sec  and  0.17  m/sec,  respectively.  If  i  -  1  ,  the  predicted  apparent 
velocltes  are  0.343  and  0.13  m/sec,  respectively.  If  l  -  2,  the  predicted 
apparent  velocites  are  0.26  and  0.14  m/sec,  respectively.  Good  agreement  is 
also  obtained  when  comparing  phases  of  t  -  2.  The  observed  phase  of  the 
differential  radius  observations  is  -82°  ±  30°;  the  predicted  values  are 
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either  0°  or  180°  for  £  -  1,  jm|  -  1;  and  ±90°  if  l  »  2,  |m|  -  2  with  an 

uncertainty  of  tl8°.  Detection  of  the  oscillation  by  an  apparent- velocity 

shift  favors  m  «•  0,  the  magnitude  of  the  observed  phase  is  only  consistent 

with  £  -  2,  and  the  sign  of  the  relative  phase  restricts  m  to  be  positive 

for  1-2.  The  160  min  oscillation  is  thus  classified  as  an  £  -  2,  m  -  2 
mode.  The  apparent-velocity  hypothesis  has  Important  implications  for 
gravity  mode  classification  (ef.  Sections  1.3.5  and  1.4.4)  and  for  the 
design  of  instruments  to  detect  ff*avity  modes. 

1 4 

1.6.3  Deviations  from  Asymptotic  Theory 

Examination  of  the  observed  differences  between  eigenfrequencies  vn+1  ^ 

and  v  . .  av  . ,  reveals  significant  departures  from  the  differences 
n,  x.  n  ,  x. 

predicted  by  asymptotic  theory.  In  particular,  the  set  of  av  n  calculated 

n  §  u 

using  the  £  -  0  eigenfrequency  spectra  of  the  5  min  modes  obtained  by 

Woodard  and  Hudson  (1983)  and  Hill  (1985a)  shows  considerable  deviation  from 

that  predicted  by  asymptotic  theory.  The  deviations  based  on  the 

observations  of  Hill  (  1985a)  from  asymptotic  theory  predictions  are 

approximately  a  factor  of  4  larger  than  any  attributable  to  errors  in 

frequency  measurements.  This  same  factor  of  4  is  also  obtained  with  the 

observations  of  Woodard  and  Hudson  (1983). 

The  scatter  in  av  „  obtained  from  Doppler  shift  observations  is  even 

n ,  t 

larger.  Based  on  the  work  of  Harvey  and  Duvall  (1984),  the  scatter  in  the 
case  for  l  >  1  is  approximately  twice  the  line  width  or  more  of  the  modes  of 
oscillation  for  n  -  20  and  four  line  widths  or  more  for  n  -  14.  These 
results  pose  what  may  be  a  very  meaningful  question:  how  can  there  be  such 
large  deviations  in  velocity  observations  relative  to  temperature 
observations  (as  inferred  from  total  irradiance  and  differential  radius 
observations)?  This  relatively  larger  scatter  appears  also  to  be  present  in 
the  velocity  observations  of  Libbrecht  and  Zirin  (1986).  The  large 
deviations  may  be  particularly  significant  in  light  of  the  reported 
discrepancies  in  rotational  splitting  that  are  obtained  using  velocity  and 
differential  radius  observations. 
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This  is  a  synopsis  of  the  work  by  Hill  and  Rosenwald  (1986). 


The  extent  to  which  the  asymptotic  theory  predictions  are  realized 
depends  upon  the  degree  to  which  the  properties  of  the  solar  interior  meet 
the  conditions  required  for  the  applicability  of  the  theory,  xs 
Chrlstensen-Dalsgaard  e£  al.  (  1  985)  state,  "The  asymptotic  analysis  is 
valid  if  the  scale  of  variation  of  the  background  state  is  much  greater  than 
all  oscillation  wavelengths  considered,  and  that  condition  is  indeed 
satisfied  throughout  most  of  the  interior  of  any  traditional  theoretical 
model  of  the  Sun.”  The  deviations  discussed  above  may  be  the  manifestation 
of  conditions  in  the  solar  interior  which  do  not  meet  the  requirements  Tor 
asymptotic  theory  to  be  applicable.  Should  this  be  the  case,  these 
deviations  take  on  a  particular  importance:  they  may  offer  a  diagnostic 
probe  to  those  regions  of  the  solar  Interior  where  a  significant  fractional 
change  occurs  in  a  characteristic  length  of  the  eigenfunction  used  in  the 
study. 

If  there  is  a  layer  in  the  Sun  that  gives  rise  to  a  sharp  change  in 
local  properties,  thus  generating  conditions  for  which  asymptotic  theory  is 
not  applicable,  the  effect  of  this  sharp  change  will  depend  on  its  location 
relative  to  the  location  of  the  nodes  of  the  eigenfunction.  The  location  of 
a  node  near  a  particular  radius  depends  on  n  and  t  and  appropriate 
simultaneous  shifts  in  n  and  l  may  keep  that  node  location  relatively  fixed. 
Thus,  if  the  model  is  a  relatively  accurate  representation  in  this  case,  it 
should  be  possible  to  demonstrate  that  the  observed  deviations  from 
asymptotic  theory  for  0  <  l  <  6  are  projections  of  a  single  distribution 
onto  seven  different  surfaces.  • 

Guided  by  this  model,  the  seven  different  sets  of  av  .  for  i  -  0,...,6 

n  ,* 

were  superimposed.  In  this  superposition  relative  shifts  were  allowed  in 

frequency,  av,  and  n,an,  which  were  proportional  to  shifts  in  t,  at.  It  was 

discovered  empirically  that  there  exists  a  au/At  and  a  an/at  which  show  that 

the  seven  apparently  unordered  distributions  of  av  .  are  projections  of  a 

n  9 1 

single  parent  function.  This  extraordinary  result  of  this  superposition  is 
shown  in  Figure  7  which  was  obtained  with 

av/al  -  -0.35;  an/ at  -  0.42  .  (12) 

The  probability  that  the  parent  function  was  obtained  from  seven  sets 
of  random  frequencies  would  appear  to  be  quite  small,  •  1.2  x  10~7.  To 
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Fig  7  The  n-dependence  of  the  superposed  spectra  of  sv  obtained  by 
the  superposition  of  the  observed  avn  for  i,  -  0,...,6  with 
shifts  in  frequency  of  av/al  -  -0.351  uHz  and  shifts  in  n  of  an/al 
-  0.^20.  The  solid  curve  is  an  estimate  of  the  quasi-periodic 
parent  function  for  that  domain  of  n  ♦  0.421  where  the  parent 
function  is  generally  overdetermined.  The  dashed  curve  is  an 
estimate  of  the  average  of  superposed  Av  .with  the  short- 
wavelength  quasi-periodic  structure  removed. 


An. 


address  the  question  are  the  quasi  per  iodi  c  results  shown  in  Figure  7 
physically  plausible,  a  numerical  experiment  was  performed.  In  this 
experiment  the  equilibrium  conditions  of  a  standard  solar  model  were  locally 
perturbed  in  a  narrow  range  in  radius  and  the  effect  of  the  localized 
perturbation  on  the  eigenfrequencies  calculated  numerically.  Excellent 
agreement  was  foind  between  the  observations  and  the  predictions  of  the 
numerical  experiment.  This  agreement  between  observation  and  theory 
indicates  that  the  type  of  mechanism  leading  to  the  observations  is 
understood  and  that  it  is  appropriate  to  use  the  results  of  the  theoretical 
analysis  to  infer  the  location  of  the  non-applicability  of  asymptotic 
theory.  It  is  inferred  that  the  effective  radius  of  this  location  is 

r/R  -  0.757  t  0.002  .  03) 

The  closeness  of  this  value  to  what  is  currently  thought  to  be  the  base  of 
the  convection  zone  suggests  that  the  norv-appl icabil ity  of  the  asymptotic 
theory  is  connected  with  the  transition  region  near  the  inner  boundary  of 
the  convection  zone. 

In  summary,  evidence  has  been  obtained  which  indicates  that  another 
important  diagnostic  tool  is  available  for  work  in  solar  seismology.  This 
probe  is  sensitive  to  those  interior  regions  of  the  sun  where  the  conditions 
for  the  applicability  of  asymptotic  theory  for  low-degree,  5  min 
oscillations  are  not  met.  Should  these  findings  be  confirmed  in  future 
work,  studies  of  the  type  described  here  may  examine  the  temporal  stability 
as  well  as  the  spatial  properties  of  the  equillbriun  conditions  of  these 
regions  where  the  scale  of  the  change  in  the  equillbriun  state  of  the  Sun  is 
of  the  order  of  the  wavelength  of  the  oscillations  used  in  the  study. 

t .6.4  Implications  of  Phas»- locked  Gravity  Modes!5 

The  work  of  Hill  U  986c)  presents  evidence  of  phase  locking  for  four 
Intermediate  degree  gravity  modes  (cf.  Section  1.4.5).  The  evidence  of 
phase  locking  of  a  single  set  of  S'avity  modes  leads  to  a  modification  of 


15.  This  is  a  synopsis  of  Hill  (1986c). 
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the  neutrino  production  rates  predicted  by  a  standard  solar  model  with  no 
solar  oscillations  present.  The  evidence  of  two  or  more  different  sets  of 
phase- locked  g*avity  modes  should  lead  to  a  periodic  modulation  in  time  of 
the  neutrino  production  rates.  The  periods  of  modulation  are  determined  by 
the  sets  of  properties  of  the  phase-locked  gravity  modes,  while  the 
amplitudes  of  modulation  are  determined  by  both  the  properties  of  the  modes 
and  the  details  of  the  neutrino  production.  Two  facts  —  that  a  modulation 
is  a  natural  consequence  of  phase  locking  and  that  the  periods  can  be 
derived  from  solar  seismological  studies  --  place  the  mode-locking 
hypothesis  in  a  unique  position  relative  to  the  many  other  suggestions  that 
have  been  offered  for  resolving  the  solar  neutrino  paradox:  a  periodic 
modulation  of  the  neutrino  production  rates  is  not  a  natural  consequence  of 
most  if  not  all  of  these  other  suggestions,  and  the  relevance  of  the  phase¬ 
locking  hypothesis  can  be  tested..  Such  a  test  can  be  made  by  looking  for  a 
correlation  between  the  predicted  periodic  behavior  of  the  neutrino 
production  rates  and  their  observed  fluctuations. 

The  presence  of  the  relatively  large  temperature  eigenfunctions  of 
phase- locked  gravity  modes  for  0.05  <  r/R  <  0.45  and  the  associated  energy 
transport  of  the  modes  are  both  expected  to  alter  the  neutrino  production 
rates.  These  two  mechanisms  lead  to  a  modification  of  the  effective 
temperature  profile  defined  for  a  given  process  in  the  Sun.  It  is 
anticipated  in  particular  that  the  former  will  be  the  more  important  process 
for  the  modification  of  the  production  rate  of  the  neutrino  from  the  decay 
of  8B. 

For  a  set  of  phase- locked  modes  belonging  to  the  multiplet  g  ,  there 

n  1 1 

exists  a  rotating  coordinate  system  with  a  rotational  speed  Q  -  *2iv' 

n  1 1 

where  each  mode  of  the  set  has  the  same  frequency  of  oscillation.  In  this 

rotating  coordinate  system,  there  will  be  one  or  more  regions  where  the 

temperatire  eigenfunctions  for  the  members  of  the  set  of  modes  are  all  in 

phase.  These  regions  of  constructive  interference  thus  all  rotate  at 

-2*v'  the  same  rate  as  the  rotating  coordinate  system.  For  a  second  set 
n,i 

of  phase- locked  modes,  an  identical  set  of  properties  exists  except  for  one: 
the  respective  rotating  coordinate  systems  will  in  general  have  slightly 
different  rotation  rates  because  the  respective  v*  t  are  slightly  different. 

It  is  the  slightly  different  rotation  rates  of  the  respective  regions 
of  constructive  interference  that  are  of  interest  here.  Periodically,  the 
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regions  of  constructive  interference  belonging  to  two  different  sets  of 
pnase- locked  nodes  will  pass  through  each  other.  Because  of  the  Importance 
of  the  higher-order  processes  indicated  by  phase  locking,  this  passage  of 
two  different  phase-locked  regions  of  oscillations  will  lead  to  enhanced 
neutrino  production  rates.  The  periodicity  of  the  enhanced  neutrino 
production  is  obtained  directly  from  the  observed  v'  (  for  the  phase- locked 
nodes.  The  results  of  the  analysis  indicate  a  typical  period  for  this 
periodicity  of  -  2  y. 

Although  the  analysis  of  the  1979  observations  for  the  properties  of 
the  interned! ate- degree  gravity  nodes  is  not  yet  completed,  the  statistical 
significance  of  the  evidence  for  the  detection  of  phase-locked  modes  is 
sufficiently  high  to  warrant  certain  prelininary  observations.  First,  the 
presence  of  the  phase-locked  gravity  modes  may  be  Important  to  the 
resolution  of  the  solar  neutrino  paradox.  More  1‘mportantly,  should  the 
node- locking  hypothesis  be  correct  in  part ,  it  can  be  subjected  to  a  direct 
observational  test,  an  opportunity  which  places  this  hypothesis  in  a  unique 
position  with  respect  to  the  other  suggestions  for  the  resolution  of  the 
paradox.  This  test,  based  on  the  observed  presence  of  several  sets  of 
phase- locked  gravity  modes  which  should  lead  to  a  periodicity  in  the 
neutrino  production  rates,  will  consist  of  looking  for  predicted 
periodicities  in,  for  example,  the  neutrino  capture  rates  recorded  at  the 
Homestake  Gold  Mine  (Rowley,  Cleveland,  and  Davis  1985).  The  predicted 
periodicities  would  be  based  on  the  1979  solar  seismological  observations. 

It  should  be  noted  that  a  positive  result  from  the  test  described  above 
would  have  an  impact  on  several  areas  of  research.  First,  of  course,  such  a 
result  would  confirm  the  presence  of  phase- locked  cavity  modes.  Second,  a 
positive  result  would  provide  the  opportunity  to  infer  in  a  quantitative 
manner  the  net  contribution  of  the  phase-locked  modes  to  the  solar  neutrino 
paradox  relative  to  other  possible  mechanisms  such  as  neutrino  oscillations. 
Third,  and  most  important,  the  confirmation  of  the  mode- locking  hypothesis 
would  demonstrate  the  value  of  solar  seismology  for  future  studies  of  solar 
neutrinos  through  the  phase-sensitive  detection  of  modulation  in  the 
neutrino  production  rates. 
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1  6 

1.7  Development  of  a  New  Observational  Technique  at  SCLERA 

The  development  of  a  new  observational  technique  was  motivated  by  two 
different  interests.  First,  an  instrument  that  could  detect  the  temperature 
eigenfunction  of  gravity  modes  on  the  solar  disk  would  augment  the 
differential  radius  observations  made  at  the  extreme  limb  and  would  be  very 
valuable  in  the  study  of  gravity  modes.  Second,  the  Influx  of  discrepant 
results  as  discussed  in  Section  1.6.1  clearly  indicates  the  desire  for  the 
development  of  new  observational  techniques  that  may  help  resolve  the 
interpretational  problems  we  currently  have. 

A  technique  to  detect  solar  oscillations  through  the  Eulerian 
perturbation  of  the  radiation  intensity./',  has  been  developed  and  tested  at 
SCLERA.  In  1985  ./’was  measured  at  48  different  wavelengths  within  the 
spectral  range  of  0.5  to  1.7  microns  (cf.  Section  1.2.1).  Since/'  has  a 
wavelength  dependence  which  is  different  from  that  associated  with  changes 
in  the  transparency  of  the  earth's  atmosphere,  the  contribution  of  4'  co  the 
observed  signal  can  be  identified.  The  contributions  were  Fourier  analyzed 
and  the  resulting  power  spectrum  indicated  to  a  high  degree  of  statistical 
significance  that  global  solar  oscillations  had  been  detected  (cf  Section 
1.5.3).  These  results  indicate  a  promising  ruture  for  the  use  of  this 
technique  in  solar  seismology  studies. 

1.8  Computer  Facilities  for  Data  Analysis  and  Theoretical  Calculations 

The  decision  was  made  in  1982  to  establish  a  dedicated  computer  system 
capable  of  meeting  SCLERA'S  off-line  computer  requirements  for  the  five 
years  1983*88.  This  decision  was  based  on  projected  real  costs  of  using 
existing  |  facilities  at  the  University  and  the  recognition  that  these 
projected  costs,  if  realized,  would  constitute  too  heavy  a  burden  on  the 
program.  During  1 980-81  most  of  the  off-line  computer  work  was  done  at  the 
University  Computer  Center;  this  work  was  funded  by  a  grant  from  the 
University  of  Arizona  in  the  amount  of  $25,000.  During  1981-82,  SCLERA'S 
computer  requirements  increased  by  more  than  100$.  The  total  cost  during 


16.  This  is  a  synopsis  of  the  work  by  Oglesby  (1986). 


this  period  was  also  met  primarily  by  grants  from  the  University.  The 
computer  requirements  for  1983~88  were  projected  to  increase  linearly  with 
time  at  the  rate  experienced  in  1981  and  1982.  However,  the  University 
instituted  in  198*1  a  major  Reduction  in  its  computing  grants.  This  made  it 
necessary  to  find  an  alternative  to  using  the  University  Computer  Center,  an 
alternative  which  would  meet  our  computing  requirements  at  a  cost  that  could 
be  supported  by  the  program.  That  alternative  was  the  decision  to  establish 
a  dedicated  computer  system  capable  of  meeting  SCLERA'S  off-line  computer 
requirements  for  .a  five  year  period.  The  five  year  plan  consisted  of 
starting  with  a  system  capable  of  meeting  the  1983-84  requirements  followed 
with  periodic  upgrades  of  the  system  as  the  program  needs  developed. 

The  computer  that  was  put  on-line  (January  1983)  is  a  Prime  250-11. 
This  initial  system  was  configured  with  a  one  Mbyte  memory,  a  64  Mbyte 
cartridge  module  disk,  a  magnetic  tape,  a  medium-speed  line  printer/plotter 
console,  five  user  terminals,  and  one  remote  dial-up  telephone  modem. 

The  projected  increase  in  computer  requirements  for  1984  proved  to  be 
accurate  and  the  first  upgrade  of  the  off-line  computer  system  was  made  in 
March  of  1984.  This  upgrade  Involved  Increasing  the  memory  to  two  Mbytes 
and  adding  a  300  Mbyte  Winchester  Disk  Drive. 

The  data  analysis  has  required  extensive  use  of  computers  because  of 
the  complexity  of  the  eigenfrequency  spectrum  and  because  of  the  relatively 
small  signals  associated  with  solar  oscillations.  The  theoretical 
calculations  have  required  extensive  use  of  computers  because  the  properties 
of  solar  models  are  being  tested  at  an  unprecedented  level  of  accuracy  in 
astrophysics.  This  complexity  of  the  data  analysis  and  theoretical 
calculations  continues  to  increase  as  projected  in  1982  as  more  and  improved 
observations  are  made.  In  1985.  the  capacity  of  the  Prime  250-11  system 
saturated;  this  was  projected  to  pose  serious  limitations  on  the  potential 
of  the  1986  program.  In  March  of  1986  the  Prime  250-11  system  was  replaced 
by  a  Celerity  1200  computer.  The  configuration  and  capacity  of  the  new 
system  is  discussed  in  Section  3.  The  Celerity  1200  with  upgrading  should 
meet  the  off-line  computer  requirements  for  several  years  with  effective  low 
operation  costs. 
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